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Abstract 

The doping of Mn in Si is attracting research attentions due to the possibihty to fabricate Si-based 
diluted magnetic semiconductors. However, the low solubility of Mn in Si favors the precipitation 
of Mn ions even at non-equilibrium growth conditions. MnSii.7 nanoparticles are the common 
precipitates, which show exotic magnetic properties in comparison with the MnSii.7 bulk phase. In 
this paper we present the static and dynamic magnetic properties of MnSii.7 nanoparticles. Using 
the Preisach model, we derive the magnetic parameters, such as the magnetization of individual 
particles, the distribution of coercive fields and the inter-particle interaction field. Time-dependent 
magnetization measurements reveal aging and memory effects, qualitatively similar to those seen 
in spin glasses . 
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I. INTRODUCTION 



Diluted magnetic semiconductors (DMS) are attracting great interest because of their 
potential use for spintronics. GaAs:Mn has recently emerged as the most popular material 
for this new technology. However, Si-based DMS would be preferably used because of the 
availability of high quality, large-size, and low-cost wafers. More importantly, the fabrication 
of Si-based DMS is compatible with the mature microelectronics technique. Based on the 
Zener model, Dietl et al. |l| predicted carrier- mediated ferromagnetism if p-type conducting 
Si is doped with 5% Mn. Using density- functional theory, Wu et al. demonstrated that 
interstitial Mn can be utilized to tune the magnetic properties of Si. Experimentally, various 
groups have reported the observation of ferromagnetism in Mn doped Si . The 

reported Curie temperatures range from 200 to 400 K. However, the opinions concerning 
the origin of the observed ferromagnetism are very diverse. Using high resolution, spatially 
resolved techniques, comprehensive material characterization reveals the clustering of Mn- 
rich phases in Mn implanted Si, namely MnSii.7 flflflSQ, which is the energetically 
most favorable Mn-silicide phase . Moverover, Mn-rich phases also form during pulsed 

laser annealing following Mn ion implantation into Si [isl. However, Mn-precipitations in 
Mn implanted GaAs are successfully suppressed [l^. Nevertheless, after considering the 
formation of MnSii.7 nanoparticles exotic magnetic properties have been observed. The 
magnetization per Mn is as large as 0.21 /i^/Mn, which is much larger than that (0.012 
/ifi/Mn) of bulk MnSii.7 j^. Ko et al. suggested the existence of multi-fold contributions 
to the observed ferromagnetism according to their observation of Mn-rich and Mn-poor 
phases . It is important to understand the dynamic magnetization of an ensemble of 



ferromagnetic nanoparticles due to its influence on technological applications. Spin-glass- 
behavior has been observed in GaAs:Mn and Ge:Mn systems containing Mn-rich clusters 



ike 
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171]. Despite numerous publications on the structural and magnetic properties of MnSii.7 
embedded in Si , information of their dynamic properties is lacking. In this 

paper, we will present the static and dynamic magnetic properties of MnSii.7 nanoparticles. 
The magnetic parameters, such as the magnetization of individual particles, the distribution 
of coercive fields and the inter-particle interaction, are deduced using the Preisach model 
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20|. Time-dependent magnetization measurements reveal a spin-glass-like behavior. 
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FIG. 1: Hysteresis loops (open circles) measured in the temperature range from 5 K to 150 K. The 
solid curves (a-d) are fittings using the Preisach model. 

II. EXPERIMENTAL METHODS 

The Mn-implanted Si samples were prepared from commercially available, Czochralski 
grown single-crystal Si(OOl) wafers, which were p-type doped with a B concentration of 
1.2 X 10^^ cm~^. Mn"*" ions were implanted at an energy of 300 keV with a fluence of 1.0 x 10^^ 
cm~^, which corresponds to a peak concentration of 0.8%, with a projected range (Rp) of 
258±82 nm. The samples were held at 350 °C during implantation to avoid amorphization. 
In order to reduce channeling effects, the angle between the sample surface normal and 
the incident beam was set to 7°. After implantation, rapid thermal annealing (RTA) was 
performed at a temperature of 800 °C for 5 min in a forming gas of N2. Magnetic properties 
were analyzed using a superconducting quantum interference device (SQUID) magnetometer 
(Quantum Design MPMS). In order to obtain the static and dynamic magnetic properties, 
we measured the magnetization depending on field, temperature and time. 
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III. RESULTS AND DISCUSSION 



A. Static magnetic properties 

The structural properties of the sample after rapid thermal annealing have been reported 



m Ref. 15|. MnSii.7 nanoparticles are formed with the average diameter of 11 nm with a sphere 
like shape. We measured the magnetization of the sample by applying the field parallel or 
perpendicular to the sample surface and did not find any difference, which hints towards 
an isotropic magnetization. However, the same effect would be seen from a collection of 
anisotropic particles with randomly-oriented easy axis. The M-H data for the sample with 
the field applied in the film plane is shown in Figure [TJ The linear diamagnetic background 
of the Si wafer has been subtracted for all shown data. The magnetization is normalized 
by the number of implanted Mn ions. It can be seen that the sample exhibits a hysteretic 
behavior with the magnetic remanence being 33% of the saturation magnetization at 5 K. 
With increasing temperature, both the coercivity (He) and remanence (Mn) drop rapidly. 
At 100 K and 150 K, the coercivity and remanence are zero. This is the evidence for 
superparamagnetism. The blocking temperature lies between 50 K and 100 K. Figure [2] 
shows the remanent magnetization measured from 5 to 200 K after shutting down the field 
of 10000 Oe at 5 K. The remanence decreases rapidly with increasing temperature and falls 
to zero at around 80 K. The coercivity (inset of Figure [2]) exhibits a similar temperature 
dependence. Note that bulk MnSii.7 is reported to exhibit weak itinerant magnetism with 
an ordering temperature of 47 K and with a very low saturation moment of 0.012 /i^/Mn 



2l|, being much different from the MnSii,7 nanoparticles investigated here. In our samples, 
the largest magnetization is 0.21 /i^/Mn. Independently, Yabuuchi et al. [8| used slightly 
different implantation and annealing conditions and realized a maximum magnetization 
of 0.2 /i^/Mn. Both values are much larger than for a bulk MnSii.7 crystal. Using first- 
principles calculation, Yabuuchi et al. have clarified that the stoichiometry, strain and charge 
accumulation as well as the interface between MnSii.7 and Si strongly influence the magnetic 
properties of MnSii.7 nanoparticles [3]. These effects well account for the experimental 
observations. 

The temperature dependent magnetization under zero field cooling and field cooling was 
measured in the following way. The sample was cooled in zero field from above room temper- 
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FIG. 2: Measured temperature dependent remanence and coercivity (inset). Both drop to zero at 
the temperature around 80 K. The red crosses display the fraction of the ferromagnetic components 
obtained by fitting using the Preisach model. Both the remanence and the fraction of ferromagnetic 
component have the same temperature dependence. 

ature to 5 K. Then a field was applied and the zero-field cooled (ZFC) magnetization curve 
was measured with increasing temperature from 5 to 300 K, after which the field-cooled (FC) 
magnetization curve was measured in the same field from 300 to 5 K with decreasing tem- 
perature. Figure [3] shows the ZFC/FC magnetization curves measured at different fields. An 
irreversible behavior is observed in ZFC/FC curves. Such an irreversibility originates from 
the anisotropy barrier blocking of the magnetization orientation in the nanoparticles cooled 
under zero field. The magnetization direction of the nanoparticles is frozen in its initial 
state at high temperature, i.e., randomly oriented. At low temperature (5 K in our case), a 
magnetic field is applied. Some small nanoparticles with small magnetic anisotropy energy 
flip along the field direction, while the larger ones do not. With increasing temperature, the 
thermal activation energy together with the field also flips the larger particles. This process 
results in the increase of the magnetization in the ZFC curve with temperature. The peak 
of the ZFC curves, Tp, is considered as the average blocking temperature of the sample. A 
notable feature is the increase of Tp with increasing applied field. Such a field dependence of 
Tp has been observed in several magnetic nanoparticle systems, such as Fe304 [22], 7-Fe203 
231], ferritin 24] and FePt 25|]. The physical origin of this behavior is discussed concerning 



the size distribution of nanoparticles 



251 ] and the inter-particle interaction 2J] 



We analyzed the experimental curves using simulations based on the Preisach approach 
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FIG. 3: Field-cooled (FC, upper branches) and zero-field-cooled (ZFC, lower branches) M-T curves 
measured at different fields. The solid lines are fitting lines using the Preisach model. 
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FIG. 4: Fitting of the hysteresis loop measured at 20 K by the Preisach model. The superparam- 
agnetic (SP) and ferromagnetic (FM) components are shown in the figure separately. 



19l . l26l | and derived the key magnetic parameters of MnSii.7 nanoparticles, namely the 
magnetization of individual particles and the distributions of coercive and inter-particle 
interaction fields. The temperature dependence of the parameters "p describing the magnetic 
properties of the nanoparticles is usually expressed by the critical temperature Tc and the 
critical exponent 7 

p = Po(l-T/Tc)^ (1) 

with p (The subscript denotes the parameters at K) substituted by the mean magnetic 
moment of an individual cluster (where 7 = F), the mean coercivity Ec (7 = Tc), or 
the dispersions Oi (7 = Fj) and Oc (7 = F'^) of the interparticle interaction and coercive 
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TABLE I: Best fit Preisach parameters. 
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field, respectively. Long-range interaction field was equal to zero. The fitting parameters 
are shown in Tab. [H 

Figure H] shows a representative fitting of the hysteresis measured at 20 K. Two com- 
ponents contribute to the measured magnetization: a superparamagnetic (SP) component 
and a ferromagnetic (FM) component with a large coercive field of 1100 Oe. Note that the 
modelled coercivity of the FM component is larger than that shown in the inset of Figure 
El and the modelled and the measured remanence are the same. The smaller coercivity seen 
in the measured hysteresis is due to the superposition of the FM and SP components. The 
fraction of the ferromagnetic component rapidly decreases with increasing temperature, and 
is plotted in Figure [21 It is worth to note that using the same set of fitting parameters we 
are able to fit the experimental curves measured at different temperatures (see Figure [ID. 

The average magnetic moment of individual particles /iq (at K) computed from the 
fitting is 4.0x10^^^ emu (around 4300 fis)- The average diameter of MnSii.7 particles in 
this sample is around 11 nm 5(]. In one particle there are approx. 21000 Mn atoms, and we 
obtain 0.204 yU^/Mn. This value is in a good agreement with the experiment (0.21 /i^/Mn). 



In the work by Yabuuchi et al. js], the samples with the same Mn fiuence were annealed 
at different temperatures to tune the size of MnSii.r particles. The ferromagnetism strongly 
depends on the average particle diameters and is maximized after annealing at 750 °C. This 
annealing temperature resulted in an average diameter of 10 nm for MnSii.7 particles. By 
the two experimental works [sl, Q] and our fitting using the Preisach approach, we confirm 
that the ferromagnetism in MnSii.y particles strongly depends on their diameters: i.e., the 
exotic ferromagnetism in MnSii.y particles is mainly a size effect [3]. 

The ferromagnetic component has a large coercive field of 1500 Oe at K, while the 
inter-particle interaction is weak. Note that ai is around 200 Oe and much smaller than the 
coercive field (Tab. [I]). These features are important for practical applications. Using the 
same fitting parameters, the ZFC/FC curves are also well fitted as shown in Figure [3l Im- 
portantly, the fitting reproduces the field dependence of Tp in the ZFC curves well. Indeed, 
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Song et al. have discussed the influence of the inter-particle interaction using the Preisach 



model 27|. In the limit of weak interactions Tp increases monotonically with applied fields, 
while in the limit of strong interactions Tp decreases monotonically with applied fields. 
Therefore, we attribute the increase of Tp at larger fields to the weak interaction between 
MnSii.7 particles. 

B. Dynamic magnetic properties 

To obtain the dynamic magnetic properties of the system, we performed time-dependent 
measurements. The time- dependent thermoremanent magnetization (TRM) was measured 
below the blocking temperature. TRM is measured by cooling the sample in an applied 
field from an initial temperature above any spin glass transition to some final temperature, 
decreasing the field to zero and observing the decaying remanent magnetization {Mr(t)). In 
our case, TRM data were taken after cooling from 300 K to 5 K and 20 K, respectively, in 
an applied field of 100 Oe. 

For superparamagnetic nanoparticles, the magnetization decay is usually exponential to 



first order, such that one should observe a single relaxation rate approximation 



28| 



Mrit) = Moe"*/" + C, (2) 

where Mq is the initial magnetization, r is the relaxation time and C is a constant. If 
the superparamagnetic nanoparticles undergo collective behavior due to direct dipole-dipo 



interaction or a large particle-size distribution, a stretched exponential form is expected |l6 |. 



Mr{t) = Moe-(*/")' + C, (3) 

where b affects the relaxation rate of the glassy component. Figure shows the TRM time- 
decays at 5 and 20 K. The stretched exponential relaxation fits better to the experimental 
data. The fitted parameters of r and b are in the typical range of a spin-glass system [ig! ]. 

In order to further confirm the glass behavior of the system, we also performed history- 
dependent magnetic memory measurements using the cooling and heating protocol suggested 



by Sun et al. 29|. We cooled the sample at 100 Oe and recorded the magnetization during 
cooling, but temporarily stopped at T = 50 K and 30 K for a period of 2 hours. During 
the waiting period, the field was switched off. After the stop, the 100 Oe field was reapplied 
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FIG. 5: TRM time-decays of the field-cooled magnetization (H = 100 Oe) at (a) 5 K and (b) 20 
K. Scattered symbols are experimental data. Solid lines (blue) are stretched-exponential fits with 
parameters of r=2120 s and 6=0.44 at 5 K, and t=1806 s and 6=0.53 at 20 K. Dashed lines (red) 
are first-order exponential fits with parameters of r=1721.2 s at 5 K and r=1720.5 s at 20 K. 





2.5 




1 1 1 1 ■ 

Memory effect 


E 

















CO 

O 


2.0 


stopj 




























1.5 




stop 1 ^v^^^^ 



1 I 1 I 1 I I 

20 30 40 50 60 

Temperature (K) 



FIG. 6: Temperature dependent memory effect in the dc magnetization. The solid squares are 
measured during cooling in a field of 100 Oe at the same rate but with a stop of 2 hours at 50 K 
and 30 K, respectively. The field is cut off during stop. The solid line is measured with continuous 
heating at the same rate after the previous cooling protocol. 
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and cooling and measuring were resumed. The temporary stop resulted in a steplike M(T) 
curve (solid squares in Figure [6]). After reaching the lowest temperature 5 K, the sample 
was heated back in the same field, and the magnetization was recorded again. The M(T) 
curve during this heating also has a steplike behavior at the stop temperature of 30 and 50 
K, then recovers the previous M(T) curve measured during cooling. The system remembers 
its thermal history. 



The observed memory effect as well as the relaxation phenomena show qua! 



ilarities to the spin glass behavior. Two explanations have been suggested 30|]. The first 



itative sim- 



one is a broad distribution of blocking temperatures originating from the distribution of the 
anisotropy energy barriers. Another explanation is the strong dipolar interaction between 
nanoparticles, which frustrates the nanomagnetic moments, and slows down their relaxation. 
Our observations support the first model. First of all, using the Preisach model we derived 
a small interaction between MnSii.7 nanoparticles. Second, the size of MnSii.7 nanoparti- 
cles approximately amounts to 11±5 nm according to the TEM observation as shown in 
Ref. la. The spin fiip time for magnetic particles depends exponentially on the particle 
size. Therefore, even a small distribution of the particle size could give a broad range of 
relaxation times. Therefore, the observed spin-glass-like behavior is attributed to the broad 
distribution of particle sizes, i.e. of anisotropy energy barriers. 



IV. CONCLUSIONS 



Static and dynamic magnetic properties of Mn implanted Si were investigated. The mag- 
netic properties of Mn-silicide nanoparticles can be well explained using the Preisach model. 
We found that there are two components (superparamagnetic and ferromagnetic fractions) 
contributing to the magnetism. The fraction of the ferromagnetic component decreases with 
increasing measurement temperature. The interaction between the magnetic nanoparticles 
is weak. Therefore, the superparamagnetism blocking-temperature increases monotonically 
with applied field. Time-dependent measurements, i.e. relaxation and memory effect, sup- 
port a spin-glass-like behavior in the investigated material system, which results from the 
size distribution of MnSii.7 nanoparticles. 
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